Excitation gaps in fractional quantum Hall states: An exact diagonalization study 
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We compute energy gaps for spin-polarized fractional quantum Hall states in the lowest Landau level at filling 
fractions V = ^, |, y and g using exact diagonalization of systems with up to 16 particles and extrapolation 
to the infinite system-size limit. The gaps calculated for a pure Coulomb interaction and ignoring finite width 
effects, disorder and LL mixing agree well with the predictions of composite fermion theory provided the log- 
arithmic corrections to the effective mass are included. This is in contrast with previous estimates, which, as 
we show, overestimated the gaps at v =2/5 and 3/7 by around 15%. We also study the reduction of the gaps as 
a result of the non-zero width of the 2D layer. We show that these effects are accurately accounted for using 
either Gaussian or 'zx Gaussian' (zG) trial wavef unctions, which we show are significantly better variational 
wavefunctions than the Fang-Howard wavefunction. The Gaussian and zG wavefunctions give Haldane pseu- 
dopotential parameters which agree with those of self-consistent LDA calculations to better than ± 0.2%. For 
quantum well parameters typical of heterostructure samples, we find gap reductions of around 20%. The exper- 
imental gaps, after accounting heuristically for disorder, are still around 40% smaller than the computed gaps. 
However, for the case of tetracene layers in metal-insulator-semiconductor (MIS) devices we find that the mea- 
sured activation gaps are close to those we compute. We discuss possible reasons why the difference between 
computed and measured activation gaps is larger in GaAs heterostructures than MIS devices. Finally, we present 
new calculations using systems with up to 18 electrons of the gap at v = | including width corrections. 



I. INTRODUCTION 

Our understanding of the FQHE [|l]] is primarily based on 
the Laughlin wavefunction (wf) and its appropriate hierar- 
chical generalizations [|[ ^, ||] for the so-called higher order 
"daughter" fractions which are many-electron wavefunctions 
in the lowest Landau level with no adjustable parameters. The 
fundamental property underlying the FQHE phenomenon is 
the existence, at certain filling fractions of the lowest Lan- 
dau level, of an incompressible ground state and an energy 
gap A in the many-body excitation spectrum. This gap is 
produced entirely by the electron-electron interaction while 
the coiTesponding non-interacting single particle energy lev- 
els are all degenerate at the particular fractional filling (i.e. 
all non-interacting single particle levels have energy /z (0^/2 in 
the lowest Landau level, where (0^ — eB/ [mc) is the cyclotron 
frequency in the magnetic field B). 

The excitation gap A is the key measure of the robustness 
of the FQHE - the incompressibility cannot be destroyed by 
weak disorder in the system if the gap is large. The behav- 
ior of the gap as a function of filling fraction in the main se- 
quence of FQHE states can also be compared to predictions 
of the composite fermion (CF) picture and used to extract the 
CF effective mass. The excitation gap at v = 1/3 has been 
theoretically estimated on the basis of exact diagonalization 
studies ^ and Monte Carlo calculations ^ P|as have the 
gaps at filling fractions v = 2/5 and v = 3/1 [|lO[. The nu- 
merically computed estimates of the gap are, however, sig- 
nificantly larger (by a factor of 2 to 3 at v = 1 /3 for exam- 
ple) than the measured gaps, A„, deduced from the activated 
temperature-dependence of the longitudinal resistivity mini- 
mum for each fraction |l2[ [l4[ . 

Here we report the results of extensive finite-size studies of 
the gap for spin-polarized excitations of electrons confined to 



the lowest Landau level (LLL) at filling fractions V = ^, |, 
I and ^ as well as detailed results for the quantum Hall state 
at V = |. We give a detailed analysis of the finite size cor- 
rections and show that previous estimates of the gap [ |lo| ] for 
the pure Coulomb interaction at v = 2/5 and 3/7 were around 
20% too high as a result of inaccurate extrapolation methods. 
For the simple exactly solvable case of a single hole in a filled 
polarized lowest Landau level, we demonstrate how an op- 
timized extrapolation scheme dramatically reduces errors in 
the estimate for the infinite system result. We compare our 
new, more accurate, results with the predictions of composite 
fermion theory [ p3[ |6, 17 1. We find that, whereas previous 
estimates were consistent with a CF effective mass which was 
independent of filling fraction, the new estimates are in better 
agreement with the CF theory which predicts a logarithmi- 
cally divergent effective mass as a function of filling factor 
as V = 1/2 is approached [|l5| |6[. The new results are also 
closer to the estimates of the effective mass from another type 
of finite-size calculations at v = 1/2 [Q. 

Other previous larger estimates of the pure Coulomb gap 
JI^ may also involve an inaccurate extrapolation to the in- 
finite system limit, but as these results were obtained using 
CF trial wf's and Monte Carlo techniques we cannot say for 
certain where the origin of this difference lies. However, we 
mention that our calculated excitation gaps are lower by as 
much as 30% than those in [19|, and some discrepancy ex- 
ists even for the pure Coulomb interaction results at v = 1/3 
where our extrapolation to the thermodynamic limit is most 
reliable. 

We clarify to what extent the discrepancy between numeri- 
cally computed gaps and those extracted from transport mea- 
surements can be attributed to finite- width effects. The large 
disagreement between experimental activation gaps Aa and 
the numerically computed gaps A^ has been an outstanding 
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problem in the subject since the first accurate measurement 
of activation gaps was reported more than fifteen years ago 
[[T4|]. There have been several previous theoretical attempts 
to compute realistic estimates of the energy gap and to iden- 
tify the source of the large discrepancy between and 
[|^, U H H H- These took account of the finite 
thickness correction (i.e. relaxing the pure 1/r Coulomb in- 
teraction approximation by including the softening introduced 
by the transverse width of the 2D layer), and of the Landau 
level mixing corrections ^7[ | . There have also been stud- 
ies of the spin-reversed excitations which are the lowest ly- 
ing excitations for small g-factors and small magnetic fields 



There are reports in the literature that the finite-width 
effects account for all the difference between measured and 
theoretically predicted gaps. Our results are at variance with 
this conclusion p2[ ] and consequently also with the results of 
p3| ] which were based on incorrect results from The er- 
ror in was originally corrected in ^5| ] . We find on the 
basis of the largest finite-size diagonalizations to date and of 
a careful analysis of the finite-size corrections that the finite- 
width corrections account for at most half of the difference 
between the computed gaps and those observed in GaAs het- 
erostructures. On the other hand the gaps observed recently in 
tetracene in metal insulator semiconductor structures are only 
slightly smaller than our estimates. We discuss the possible 
reasons for these discrepancies. We argue that they are un- 
likely to be due to spin-reversed excitations or Landau-level 
mixing and suggest that they are the result of disorder effects 
which may affect the activation energy for transport differ- 
ently in heterostructures and MIS devices. 

We show that it is unlikely for there to be a transition from 
an incompressible to a compressible state at fixed filling fac- 
tor, for V = 5' fi 7' caused by a gap collapse induced entirely 
by the softening of the Coulomb interaction due to the finite 
thickness corrections. Such a transition has been conjectured 
to occur in the second Landau level [ |3l| , where the FQHE 
is much less robust. It may also happen in situations where 
increasing the width in the transverse direction changes the 
symmetry of the sub-band wf [^. Alternatively, a new kind 
of FQH state can arise in square, parabolic or double wells, 
where, for large enough well width, the wf may split into an 
effective double layer structure at the two ends of the well 
with a central self-consistent barrier separating these two ef- 
fective layers [M]. In the regular GaAs heterostructure system 
[|ri|, |l^ |ljwe find the lowest Landau level v 5, §, ^ 
FQHE to be robust with respect to the finite thickness effect, 
with A > even for the largest possible (and physically al- 
lowed) transverse thickness. However the actual value of A 
may become rather small and one might have to go to very 
low temperatures (and very high quality, low disorder sam- 
ples) to observe the FQHE. Our results are in conflict with the 
claim by Park et al. 
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35|] that the finite width alone 
can lead to the loss of incompressibility at a filling fraction 
V = p/{2p + l) for some finite value p — pc- 

We also compare the various approximate methods for 
accounting for finite thickness effects based on interfacial 
trial wf's with those taken from self-consistent local den- 



sity approximation (SCLDA) calculations P4 , pq ]. Previous 
model calculations have used Gaussian and Fang-Howard en- 
velope wf's and the Zhang-DasSarma (ZDS) model interac- 
tion |20|. We introduce a new variational envelope wf, the 
"z X Gaussian" (zG). We find that both the zG and the Gaus- 
sian envelope wf's give Haldane pseudopotential parameters 
which agree to within fractions of a percent with those from 
the full SCLDA wf's with the zG giving slightly more accu- 
rate results at the densities used in experiment. However, both 
give essentially indistinguishable results for excitation ener- 
gies and gaps from those taken from the SCLDA wf's. This 
result shows that accurate finite-size studies of finite-width ef- 
fects require only the determination of the appropriate width 
parameter in either the Gaussian or zG description and do not 
require the use of SCLDA based tables of parameters [^2|]. 

We show that depending on the sub-band density either the 
Gaussian or zG variational wf provide substantial quantita- 
tive improvements over the well-known Fang-Howard varia- 
tional wf [|7]| which has been employed extensively in het- 
erostructure electronic calculations. Indeed, it turns out that 
the Fang-Howard wf generally overestimates the kinetic en- 
ergy, and consequently predicts significantly too large width. 
The expectation value of the energy and other quantities of in- 
terest in this context can be calculated analytically for these 
variational wf's, and in the case of the Gaussian, it is easy to 
perform expansions for either very small or very large width. 
Finally in the Appendix, we explain why the ZDS model is 
not reliable directly for predicting finite thickness corrections, 
but we present a simple modification which corrects its main 
shortcoming. 

The remainder of this paper is organized as follows: In 
section II we describe the diagonalization of the A^-particle 
Hamiltonian in the spherical geometry and give the definitions 
of the quasiparticle, quasihole and gap energies. In section III, 
we discuss the extrapolation to the N ^ °° limit and in section 
IV we compare the variation of the calculated gaps with fill- 
ing fraction v with the predictions of composite fermion the- 
ory [|l5} |l7|]. In section V we show how variational wf's 
can be used to model finite width effects. In section VI, we 
compute the reduction of the energy gaps as a function of the 
finite width and in section VII we compare our results for the 
gap energies with experimentally reported estimates of gaps. 



II. QUASIPARTICLE AND QUASIHOLE ENERGIES 

We model the two-dimensional electron gas using Hal- 
dane's spherical geometry [||]. Particles with coordinates 
(7?,9i,(|)i) move in a monopolar magnetic field of strength 
B = SH/elf' which gives rise to 25+ 1 linearly independent 
cyclotron orbits in the lowest Landau level. The single parti- 
cle orbitals on the surface of the sphere for the particles in the 
lowest Landau level are then functions x|/(9, ,(|),) which are the 
lowest energy eigen states of the kinetic energy. 

In the lowest Landau level, the interaction between particles 
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is written 

v(/7)=LLk,^™(o-) (1) 

m i<j 



where Pm{ij) projects onto states in which particles / and j 
have relative angular momentum nil and y„, gives their inter- 
action energy for this relative angular momentum. The set Vm, 
called Haldane pseudopotentials [|]|, completely characterizes 
the interaction between particles confined to the lowest Lan- 
dau level. In terms of the electron-electron interaction, V{r), 
they are defined in the plane by [^sj 

(2) 

where n refers to the Landau level and V{r) stands for the 
electron electron interaction. The corresponding integrals for 
electrons on the surface of a sphere are described in [^]. In 
the lowest Landau level n ~Q, the first Laguerre polynomial 
in equation (Q) is equal to unity. As we shall discuss in the 
next section, the effect of the finite width of the wf, (j) is 

incorporated in these pseudopotential parameters Vm \ In the 
following, we will drop the superscript n = and denote the 
Haldane pseudopotentials for the lowest Landau level by V,,,. 

The method for computing excitation energies and gaps in 
thisgeometry has been described in detail in many places [|[ 
0, |o[. According to the hierarchy model, the FQHE ground 
states at filling fraction v occur for a system of particles 
when the total flux 25 is given by 
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energy of hole at p=^ 
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Figure 1: The energy of a single hole. The asterisks show e^{N) and 
the circles show the corrected energy e^[N). The extrapolation in 
1/A' of the corrected energies pass through the exact value 
whereas extrapolation of £^{N) in 1/A' would give incorrect results. 

while at filling fractions v ~ p / {2p + I) with p 7^ 1 we have 
to estimate E2Sq (N) by interpolating (or extrapolating for the 
largest system sizes) between system sizes for which we can 

compute £250 (■^) ■ 



2SQiv,N)^v-^N + X{v) (3) 

where X{v) is the shift function Jiofa, which is a character- 
istic of the geometry of the system (in this case the sphere) 
[^]. Laughlin's [Q] elementary fractionally charged exci- 
tations from the FQHE ground state at filling fraction v = 
p/{2p + 1) coiTespond to the ground state configuration of 
a system with additional/missing flux ± 1 /p, 

25±i/p(v,A^)=25o(v,A^)±-. (4) 

p 

At V = 1 /m there are systems with both 250 and 2S±i/p both 
integer for all integer A^. At other filling fractions 2Sq and 
25-1- i/p are never both integer for the same A^. For example 
at V = 3/7, 25o is integer when the particle number is A^ = 3n 
(n integer) while 25-|-i/3 is integerforA^ = 3«=F 1, respectively. 
We take the energy to nucleate a single quasiparticle/quasihole 
in a system of N particles at filling fraction v, (N), to be the 
total energy difference between the lowest energy state with 
total flux 2S±i/p{v,N) and the total ground state energy the 
system would have at 2So{v,N) for the same A^, i.e. 

e^{N)^E2s^,/„{N)-E2SoiN). (5) 

Here E2s±i/p{N) is the total energy of the system of A^ par- 
ticles in their ground state in 2S±i/p flux quanta. For fill- 
ing fractions v = 1/m we can calculate both energies directly. 



III. ENERGY GAPS: FINITE SIZE EFFECTS AND THE 
THERMODYNAMIC LIMIT 

From studying the variation with system size of the ener- 
gies to nucleate quasiparticles and quasiholes in finite-size 
systems we estimate the excitation energies in the thermody- 
namic limit. It is essential that the extrapolation procedure is 
carried out carefully. Firstly, if working with the Coulomb in- 
teraction it is usual to quote energies in units of e'^/ei'^, where 

£q = ^Jhc I eB is the magnetic length and e is the dielectric 
constant for the medium. Howev er, for systems with n umber 
density on the sphere, £q = -^Z ( 1 / (27tn5) ) (A^ / (25) ) and so 
for systems at fixed density, the magnetic length depends on 
the particle number and the total flux through the ratio N /2S. 
In order to compare quantities measured in the same units we 
convert all energies by using the magnetic length in the infinite 
system = ^Jv / {2Tins) . 

There is also a systematic contribution to the excitation en- 
ergy in a finite size system which scales to zero in the thermo- 
dynamic limit, which we can take account of explicitly JlO|]. 
When the localized quasiparticle/quasihole excitation which 
is formed in a system of A^ particles around the point on the 
sphere RQ. with £2 a unit vector pointing away from the ori- 
gin a charge ±qe with q= l/(2/:> + 1) is concentrated around 
this point. This charge has come from the rest of the system. 
There is then a contribution, Aq, to the energy of the system 
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from the non-uniform distribution of charge on the surface of 
the sphere which, in units of e^/e^o, is given by 



'1 



(6) 



To extrapolate to the infinite system size Hmit it is better to 
remove this contribution expUcitly and study the corrected 
quasihole and quasiparticle energies 



^{N)^e^{N)-A,{v). 



(7) 



We also define the corrected gap energies to be the sum of 
quasiparticle and quasihole energies 



(8) 



We denote the limit N ^ °° of the gap and quasihole, quasi- 
particle excitation energies by e|f ' and respectively. 

To illustrate the importance of working with these corrected 
energies, we show results for a single hole at v = 1, which is a 
case for which the energy can be computed analytically using 
the exact expression for the energy of a filled Landau level 
M\. We find 



2N+1 



1 



(9) 



The contribution — 1 / y/lN is just the correction A i ( 1 ) for the 
case V = 1. Both e'l{N) and e^{N) are shown as a function 
1 /N in Figure [l|. It is clear from the figure that extrapolation 
of e^{N) with 1 /N would give spurious results because of the 
contribution of which varies as 1/^/N. By contrast, 

extrapolation of (N) with 1 /N gives the correct result 

In Figure ^ we show If^p (N) and e^p (N) and their extrapo- 
lations toN = °° using least squares fits to linear and quadratic 
functions in 1 /N. We take the difference, 0.0005 for the quasi- 
particle and 0.0002 for the quasihole, in the estimates from the 
two different extrapolation procedures as our measures of the 
accuracy of the extrapolation. In fact, inclusion of the 1 /N^ 
term in the fit does not improve the value significantly. So, 
in the following, we will always use linear extrapolation in 
1 /N to compute gaps in the thermodynamic limit. Figure || 
shows the gap energies at v = 1 /3 and 2/5 as functions ofl/N 
and the extrapolations to N = °° limit together with the esti- 
mates based on the study of trial CF-wf's [p^. Plotted are the 
sum of quasiparticle and quasihole energies e^{N) +e^{N), 
using the correction A^(v) and the energy eexc{v) of the 
neutral excitation with L ~ L,„ax (cf. caption to Figure ||), cor- 
responding to maximum separation of the quasiparticle and 
quasihole on the sphere, again corrected by the term Ag{v) 
which stands for the Coulomb energy between the quasipar- 
ticle and quasihole. As can be seen, the size dependence of 
the exciton energies is much less smooth than that of the sum 
of qp- and qh-energies. Indeed, if they were known only for 
small systems, extrapolation to the bulk limit would be inac- 
curate. Only for the largest systems, does the size dependence 
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Figure 2: The quasiparticle (ejyj(A')) and quasihole (e^yj(A')) en- 
ergies at V = 1/3 and the extrapolation using linear and quadratic 
functions of 1/A'. We take the small differences in the extrapolated 
result as a measure of the accuracy of the extrapolation. 



of exciton energies become smooth and allow reUable extrap- 
olation to the thermodynamic limit, which is consistent with 
that based on the sum of qp- and qh-energies, although less 
accurate (see Figure ||). 

Figure ^ shows the quasiparticle and quasihole energies at 
V = 3/7 and 4/9. Although the estimate at v = 1/3 is close to 
the values quoted previously [p^, the values at v = 2/5 and 
3 /7 are around 20% smaller although still within the large un- 
certainties of the earlier calculation. Our latest estimates are 
more accurate as a consequence of a better understanding of 
finite size effects in addition to being able to diagonalize the 
Hamiltonians for larger systems (with up to around 100 mil- 
lion basis states) - 15 particles instead of 11 particles for the 
quasiparticle and quasihole at v = 2/5 and 16 particles instead 
of 13 for the quasiparticle at v = 3/7. It is interesting to note 
that the extrapolation of quasiparticle and quasihole energies 
at V = 1/3 based on small sizes (A^ = 4,5,6) yields the values 
£j/3 w 0.0757, e+/3 « 0.0267 and for the gap A1/3 w 0.1024. 
These are within about one percent of our best estimates of 



the bulk limit of = 0.0749, 6^/3 = 0.0263 and for the 
gap Ai/3 — 0.1012, obtained using system sizes up to A' = 12 
and performing the extrapolation by linear polynomial fit in 
1/A^. Likewise at v = 2/5, extrapolation using the results at 
= 5 , 7 , 9 yields values for the bulk limit for the quasiparticle 
and quasihole energies of Ejyj ~ 0.0431 and ej^^ w 0.00920 
and a value for the gap A2/5 ~ 0.0523, while our best esti- 
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Figure 3: The gap energies e^{N) at v = 1/3 and 2/5, computed 
using (small solid dots), and the neutral exciton energies for L = 
A'atv=l/3andL=(A' + 2)/2atv = 2/5 (large shaded dots) . Also 
shown by circles on the vertical axes are the estimates of the gap 
energies of Jain and Kamilla [|l9| obtained from an analysis of trial 
composite fermion wf 's. The straight lines denote the best linear 
(in 1 /N) fit to the data points. The intercepts give the estimate of 
the gap energy e^(°°) neglecting corrections due to non-zero width 
effects and higher Landau levels. 
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mate based on system sizes 7 < N < 15 are E2/5 = 0.0398, 
^2/5 ~ 0.0102 and the gap A2/5 = 0.0500, corresponding to 
a difference for the gap of about 5 percent. This observation 
makes us confident that it is now also possible to compute re- 
liable buUc limit values for the gaps at V = 3/7 and 4/9. Our 
values are A3/7 = 0.035 and A4/9 = 0.027. The systems at 
V = 4/9 were inaccessible in our earlier work [ 10|. 



IV. EFFECTIVE MASS OF COMPOSITE FERMIONS 

Our estimates of the gap energies in the sequence of states 
V = p/ {2p + 1 ) are compared in Figure [ill with the predictions 
of CF theory |l|, 0, which for this sequence gives (in 
units of e^/e£o) 



1 



2 |2p+l|(ln|2p+l|+C'; 



(10) 



Choosing C' = 4.11 to fit the gap at v = 1/3 gives the gaps at 
V = 2/5, 3/7 and 4/9 to be 0.0549, 0.0371 and 0.0276 which 
are remarkably close to the estimates we obtain. We also note 
that, whereas the earlier estimates were better fitted by assum- 
ing that the gaps were simply proportional to 1/ {2p + 1) (i.e. 
ignoring the logarithmic corrections). Figure ^ shows that the 
new results are better described by the theory when the log- 
arithmic corrections to the gap are included. This translates 



Figure 4: The quasiparticle and quasihole energies (N) at v = 
3/7 and 4/9 and the best linear fits to these points. The sum of the 
two linear functions can be taken as a measure of the gap energies 
e^{N) (these cannot be computed directly for these filling fractions 
as the systems with single quasiparticle and quasiholes have different 
numbers of particles). 



into an effective mass 



m*ip) ^n^ 



F{p) 



where 



F{p) = -[ln|2p+ 11+4.11]. 

7t 



(11) 



(12) 



The effective mass of CF's has also been estimated by 
studying the variation with system size of the ground state 
energy for systems of electrons close to v = 1 /2 on the sphere 
with 25o(l /2,N) = 2A? - 2 [|l8l|. These studies gave F - 5, 
which is about 25% larger than the value we obtain for p — A. 
One would expect that, in systems close to v = 1/2, the ef- 
fective mass would be larger than at v = 4/9 but still finite as 
the long-wavelength fluctuations of the Chern-Simons gauge 
field, which give rise to the logarithmic divergence in the ef- 
fective mass, will be cut off by the level spacing. 

In [10 1 estimates of the gap energies based on collective 
excitations were also presented. On a sphere the effective 
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The estimates based on our finite-size studies of systems with v = 
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corrections from ([lo|). The constants C and C' are chosen to give the 
correct gap at p = 1 (v = 1 /3). 



wavevector of a collective excitation is kgj^=L/R where L is 
the total angular momentum of the system. This lowest-lying 
collective excitation should correspond to a well-separated 
quasiparticle-quasihole pair in the limit of large L. In the hi- 
erarchy picture, the separation of the particle and hole should 
be 2RL/N, so the maximum separation possible occurs when 
L = Ni. Here A^j is the number of particles in the condensate 
of the highest (f th) level of the hierarchy that occurs: A'o = A^, 
Nx = {N + 2)/2, N2 = (A^ + 6)/3 and N3 =JN+ 12) /A- for 
V = 1/3, 2/5, 3/7 and 4/9 respectively [| |, |l^]. Extrap- 
olations to the infinite system limit of the L—N excitations 
should therefore give an estimate of the gap energies. The re- 
sults for V = 1/3 and 2/5 are also included in Fig ^ It is clear 
from the figure that an extrapolation based on the exciton en- 
ergies would not be as smooth as that based on the charged 
excitations. 

We believe that the exciton energies in the small systems 
accessible to direct diagonalization are not as reliable a basis 
for extracting estimates of the gaps as the sum of the quasihole 
and quasiparticle energies. The principal reason for this is that 
the quasihole and quasiparticle states are actually ground state 
configurations of particles in total flux 2S±iip and they are 
well-separated in energy from all excitations. On the other 
hand, although the neutral excitations are minimum energy 
states for the quantum numbers concerned, they are close to 
the continuum of excitations for these quantum numbers and 
this gives scope for large finite size effects, in addition to lead- 
ing to poor convergence and numerical instability of vector 
iteration (Lanczos type) diagonaUzation methods. With the 
possible exception of the systems at v = 1/3, it is also clear 
that the system sizes accessible to direct diagonalization are 
not large enough to accommodate two excitations without sig- 
nificant overlap of the charge profiles of the quasiparticle and 
quasihole. In Figure ||, we show the density profile of the 14 
particle exciton at v = 2/5, with the corresponding quasipar- 
ticle and quasihole density profiles for a 13 particle system 
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Figure 6: Density profiles of excitations at v = 2/5 as a function of 
great circle distance from the north pole: The L = 8 collective excita- 
tion (exciton) for the 14 particle system, the quasihole (L — 7/2) and 
the quasiparticle (L = 4) for the 13 particle system. The projection of 
angular momentum is maximal {L, = L) in all cases. The origin for 
the quasiparticle has been shifted, so that the point at the south pole 
coincides with that for the exciton. In the thermodynamic limit, the 
exciton with these quantum numbers becomes a quasihole localized 
at the north pole and a quasiparticle localized at the south pole. It is 
clear that, even with 14 particles at v = 2/5, there is still significant 
overlap between the density variations associated the quasiparticle 
and the quasihole localized about opposite poles. This probably ex- 
plains the large finite size effects seen (Figure |3|) in the exciton energy 
as a function of A'. 



overlaid for comparison. The quasiparticle and quasihole at 
the opposite poles are clearly visible, but the system is not 
large enough for the density profiles not to interfere. 



V. INTERFACIAL WAVEFUNCTION AND MODIFIED 
INTERACTION 

The finite width of the sub-band envelope wf in the direc- 
tion perpendicular to the plane of the two-dimensional elec- 
tron gas can be incorporated into an effective interaction be- 
tween electrons in the plane. With the magnetic field perpen- 
dicular to the plane the single-particle orbitals can be written: 



(13) 



The in plane wf's \\i(x,y) are eigenfunctions of the free elec- 
tron Hamiltonian in a perpendicular magnetic field, while ^(z) 
satisfies the Schrodinger equation for a particle in the confin- 
ing potential of the quantum well or heterostructure [ p2[ ] . The 
effective interaction between particles V ( | n — r2 1 ) at positions 
ri = {xi,yi ) and r2 in the two-dimensional electron gas is then 
given by 



V(|ri-r2|) = (eVe) J dzi J 



dZ2 



(14) 

The study of finite-size systems is based on exact diag- 
onalization or the study of variational trial wf's for parti- 
cles in a given Landau level with the interparticle interaction 
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taken to be V(|n —r2\)- For particles on a sphere, the in- 
teraction V{\ri — r2|) projected onto a given Landau level is 
characterized by Haldane's pseudopotential parameters {V,,,} 
{m = 0, 1, . . .)■ Once these are known the exact diagonaliza- 
tion proceeds exactly as in the zero-width case. (We note that 
as the pseudopotential parameters are computed from the ef- 
fective interaction, which assumes a planar geometry, there is 
no attempt to account for any effects of the curvature of the 
sphere on the finite width effects.) 

Within a local density functional scheme the wf ^(z) satis- 
fies the equation 



d 



1 



2 dz m* (z) dz 



where V^ff includes the effect of the confinement poten- 
tial (including the effect of the depletion layer), the Hartree 
self-interaction and exchange-correlation. For GaAs-GaAlAs 
quantum wells the jump in m* and the dielectric constant, e 
across the interface are small and to a good approximation 
both quantities can be taken to be independent of of z (see 
Table 1 in [^6|]) and the equation simplifies. In [ p^ this equa- 
tion was solved numerically for various geometries and the 
results presented in the form of tables of pseudopotentials for 
quantum wells and heterostructures for various values of the 
electron density and device parameters. Here we show that 
the values of the pseudopotentials characterizing the Coulomb 
interaction in the finite-width geometries can be very accu- 
rately computed using a Gaussian and a new trial wf, the 
'zx Gaussian' (zG), thereby allowing the effect on the pseu- 
dopotentials of a finite-width in the direction perpendicular to 
the 2D electron gas to be encoded in just one variational pa- 
rameter, ie the width of the wf, parametrized by the standard 
deviation w of the probability density. 

The self-consistent computation of ^(z) and Ve//(z) is stan- 
dard and follows the procedure given in The poten- 
tial Veff{z) is written 



Veffiz) = Vwiz)+VHiz)+Vxc{z) 



(16) 



where Vw (z) is the confining potential of the quantum well or 
heterostructure (including image charge effects and the effect 
of the depletion layer) and Vxc is the exchange-correlation 
potential 



Vxc = [1 +0.7734xln(l +x-^)] 



R 



(17) 



where p = (4/971)'/^ x = r,/21, r, = (47^:fl*«(^)/3)-'/^ with 
a*and R the effective Bohr radius and Rydberg in GaAs. The 
Hartree potential is given by 



semiconductor there will be rijgpi (= NaZd) charges per unit 
area of the interface distributed evenly across the depletion 
layer of width zd- 

We obtain ^(z) by solving ( [l5| ) using trial forms for ^(z) and 
cornpare the results with those obtained by numerical solution 
in [|2[. The trial waveforms we have studied are the Fang- 
Howard (FH), which is zero for negative z and for positive z is 
given by 



Vh{z) 



^jdzj rfz'|z-z'| (IC(z)l'-p(z)) (IC(Z')I'-P(^')) 

(18) 

where p(z) is the (neutralizing) charge density of the doping 
ions which are taken to be far way from the interface. In 
Vh is referred to as the potential due to the induced charges 
or Vs- In the presence of Na acceptors per unit volume in the 



C(z) ^ zexp(-/7z/2), 



the Gaussian 



C(z) exp(-(z- aw)^/4w^), 



(19) 



(20) 



and the 'zx Gaussian' (zG) wf, which is again zero for nega- 
tive z and for positive z is given by 



°= zexp(-z^/9c^) 



(21) 



The width W of these wave functions can be characterized by 
the standard deviation of the corresponding probability den- 
sity. It is given in terms of the parameters i>, w,c as follows. 



b 



for the Fang-Howard wf. 



for the Gaussian wf and 



Wg = w 



3V8 + 37t- 16 
2VS ^ 



1.01016 xc 



(22) 



(23) 



(24) 



for the zG wf. 

We determine the parameters b, w, a and c variationally. 
We have found that, expanding the expression for Vxc in 
( p^ ) in X and keeping only the constant and linear terms, 
reproduces the correct expectation value for the exchange- 
correlation energy to within 0.1% for all three trial wf 's, while 
including the quadratic term affects only the fifth significant 
figure. The small x expansion works well because the dom- 
inant contribution to the exchange-correlation energy comes 
from the region in which the density is high (x small). Us- 
ing this expansion allows us, for the three trial forms for ^(z), 
to compute analytically all the integrals involved in computing 
Veff{z) and, hence, also the expectation value for the sub-band 
energy E. For the case of the Gaussian trial wf, the effective 
interaction ( p4[ ) can be written in closed form in terms of the 
Bessel function Kq, 



Voir) 



1 



(25) 



In Figure ^ we compare the estimates (see |l5|) of the sub- 
band energy for the three variational wf's: ^(zG), £ (Gauss) 
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Figure 7: Comparison of the variational estimates for the sub-band 
energies (see llSl) and of the standard deviation of the charge dis- 
tribution (width) for the tiiree variational wf 's. The results for the 
width from direct numerical solution of the equation by Ortalano 
et al [ p^ ] are also included. The top panel shows the ratio of 
the variational estimates ii(FH)/£(zG) and the second panel shows 
£(Gauss)/ii(zG). The depletion layer density, n^^pi, is quoted as a 
fraction of the electron density in the subband. 



and ^(FH) . For higher densities (jis ^ lO'Vcm^), ^(Gauss) 
gives the lowest variational estimate, while for lower densities 
^(zG) gives the lowest estimate. For all densities in the range 
we have studied, we find that the Fang-Howard wf is worse as 
a variational wf than the zG and significantly worse at higher 
densities than the Gaussian. This is because the FH wf has 
too high a kinetic energy which it can only reduce by spread- 
ing the density wider. Although the variational estimate of the 
energy for the FH wf differs by a factor which only varies be- 
tween 5% and 20%, the width of its distribution, as measured 
by the standard deviation, is significantly larger (~ 50%) than 
for the other two wf 's. 

Given that the integrals involved in using the Gaussian or 
zG wf 's can be performed analytically and are more accurate 
as trial wf 's, it is perhaps surprising that these wf 's have not 
been more widely used in the study of heterostructures and 
quantum wells. Of the two, the Gaussian is easier to use, al- 
though it will be less well-adapted to MOS devices with large 
band gap discontinuities. For the heterostructures considered 
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Figure 8: The Haldane pseudopotentials for the interaction Vc(|Fi — 
F2I) ( ps] ) projected onto the lowest and second Landau levels as a 
function of angular momentum L for particles in heterostructures 
with densities n. The results of direct numerical solution of equa- 
tion ( |r^ ) taken from are shown as crosses and results based on 
the Gaussian trial wf for ^(z) as circles. The width parameter w of 
( p^ for the trial wavefunction is chosen so that the Haldane pseu- 
dopotential V\ matches that obtained by numerical integration. The 
differences between the results based on the Gaussian interface wave 
function and the numerically computed local density approximation 
is at the fraction of a percent level. 



below we use a (conduction) band gap discontinuity of 200 
mV - the value appropriate for a GaAs/GaAso.66Alo.33 het- 
erostructure. On the other hand, the 'zG', should become 
more favorable as a variational wf, when the band gap dis- 
continuity is large and the effect of the boundary is well ap- 
proximated by a hard wall. 

Figure ^ shows the Haldane pseudopotentials for the inter- 
action V{\r\ —r2\) projected onto the lowest and second Lan- 
dau levels for heterostructures with densities appropriate to 
samples studied experimentally. Here, we determine the width 
of the trial interfacial wf 's by requiring that the Haldane pseu- 
dopotential V\ from is exactly reproduced. We note here, 
that results for the pseudopotentials in [ p2| ] were for a value of 
the magnetic length to which coincides with the Bohr radius 
a*g sa lOnm of electrons in GaAs. The results for the second 
Landau level were for the density = 6 x lO'*' cm^^ used in 



| |22| ] . The study of the second Landau level in was moti- 
vated by the results reported in at filling fraction v = 5/2. 
However, the interfacial wave function is determined by the 
total number of electrons, which for the sample studied in 
[ |4^ ] was w 3 X IQ^^cnT^, and not by the fraction occu- 
pying the second LL («i'' = 6 x lO^'cm^^) incorrectly used 
in p^ ] . For this reason, the conclusions regarding the v = 5/2 
state of Reference are incorrect. 

It is clear that the use of the Gaussian trial wf yields results 
which are essentially indistinguishable from the results of the 
exact numerical solution for ^(z). We find very similar results 
for the zG. In Figure ^ we show the difference between the 
pseudopotentials computed exactly by solving numerically for 
the interface wf ^(z) (taken from [E2|]) and those obtained us- 
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Figure 9: Errors in the Haldane pseudopotentials computed using the 
Fang-Howard, z x Gaussian (zG) and Gaussian interface trial wf 's for 
two different densities. The comparison is with the values reported 
in p2|]. The variational parameters are determined such that Vi is 
correctly reproduced. The zG leads to errors for V,„ I which are 
roughly half as big as those with the Fang-Howard wf. The Gaussian 
wf is even better in reproducing the LDA results with a maximum 
error of less then 0.2 percent (for Vq)- The width parameters for the 
three variational wavefunctions (given above each curve) are approx- 
imately equal for both densities for all three wavefunctions, implying 
that fixing the pseudopotential Vi is effectively equivalent to fixing 
the standard deviation of the charge distribution. 



ing the Gaussian and FH wf s. The errors obtained using the 
FH wf are at the 1% level while those obtained for the Gaus- 
sian are at the 0.1% level. Those obtained using the Gaussian 
trial wf are smaller than other uncertainties in the model such 
as those related to the value chosen for the depletion density, 
ndepi- Finite-width effects on the Haldane pseudopotentials 
are clearly accurately captured by the Gaussian (and zG) trial 
wf 's. Given the fact that the pseudopotentials V,,, only depend 
on the width parameter (w for the Gaussian wf, b for the Fang- 
Howard wf and c for the zG), it is clear that the use of these 
trial wf 's massively simplifies the study of finite-width effects 
when compared to the numerical integration of ( [l5| ) and tab- 
ulation of pseudopotentials used in [p2|]. For the case of the 
Gaussian, we also have an analytic expression for the effective 
interaction V (|n — r2|), cf equation (p^). 

The tables 1-lV of reference can be summarized by list- 
ing the effective width of the Gaussian interface wf for which 
the Haldane pseudopotential of order m = 1 is exactly repro- 
duced. In Table | we list the width parameters for all tabulated 
cases. 



I 

parabolic QW 


II 

Heterointerface 


III 

square QW 


IV 

Heterointerface 


«5 w [nm] 
0.49 19.3813 
0.60 24.1478 
0.73 28.9638 
0.85 33.1890 


ris w[nm\ 
0.10 9.30690 
0.50 7.13556 
1.00 5.76001 
2.00 4.61801 
3.00 4.06015 


ris w[nm\ 
0.10 3.21628 
0.50 3.21854 
1.00 3.22610 
5.00 3.31037 


ris w[nm\ 
0.60 6.72784 



Table I: Width parameters for all tabulated results of Reference [[22|]. 
The electron density n, is measured in units of lO"/™^ while the 
width w of the interfacial wf is given in nm. 
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Figure 10: The quasihole energy e^yj(A') and quasiparticle energy 
e^^3 (A') and the best linear fits to these points computed as a function 
of the width w of the density distribution computed using Gaussian 
trial wf 's. The sum of the two linear functions can be taken as a 
measure of the gap energies ^^^^(A'). 



VI. FINITE WIDTH EFFECTS ON ENERGY GAPS 



A. Filling fractions v : 



p 



With these modified potentials we have repeated the cal- 
culations described in section 2. Using the Gaussian wf 
parametrized by its width w we compute the Haldane pseu- 
dopotentials as a function of w. By exact diagonalization just 
as in the pure Coulomb case of Section 3, we compute width 
dependent excitation energies for all possible system sizes and 
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0.12 




0.0 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

width of interface wave function w [Iq] 

Figure 1 1 : Estimates of the energy gaps in tlie tliermodynamic limit 
as a function of the width of the subband wf (taken as the standard 
deviation of the charge distribution). The solid lines show the fits to 
the interpolation formula (cf. Equation (Gq) and Table 



V 




(|) [degrees] 


a 


b 


1/3 


0.1012 


34.18 


0.1468 


1.542 


2/5 


0.0500 


36.07 


0.1935 


1.866 


3/7 


0.0350 


34.98 


0.2078 


1.851 



Table II: Interpolation function for the gap energy as a function of 
width (p6|): parameters for filling fractions v = 1/3, 2/5, 3/7. 

perform for each value of the width parameter w an extrapo- 
lation to the thermodynamic limit N ^ °°. As an example we 
show in Figure |o| the size and width dependent quasiparticle 
^^^{N) and quasihole energies e^p{N) at v = 1/3. For each 
width w we use linear extrapolation in 1 /N to estimate the gap 
energy in the thermodynamic limit as a function of width. The 
size dependence at finite width is qualitatively the same as at 
w = 0. This same procedure was also employed for the calcu- 
lation of width dependent quasiparticle and quasihole energies 
at V = 1 /3 and v = 2/5, and the corresponding energy gaps in 
the thermodynamic limit. The result of these calculation are 
shown in Figure [ll|. The full lines in Figure [ll] correspond to 
interpolation functions of the form 



Carlo simulations (MC) of CF trial wavefunctions, which as 
mentioned in section TV, give larger gaps than our results for 
the bare Coulomb interaction even at v = 1/3. This discrep- 
ancy exists throughout the range of w/£o in the figure with 
our estimates being between ^5% smaller (for v = 1 /3) and 
~ 25% smaller (for v = 3/7 and 4/9). (It is not surprising 
that the difference does not depend strongly on w/Iq: While 
the energies are affected by the width w through the varia- 
tion of the effective interaction, the wavefunctions are not ex- 
pected to change significantly [ p[.) The gaps as a function of 
width have also been estimated [ [43| ] using a field theoretic ap- 
proach [^, which constructs explicit CF wavefunctions out 
of Chern-Simons gauge-transformed fermions. Energies of 
ground and excited states can be computed within this theory 
at the Hartree-Fock level. The theory needs to cut off the in- 
teraction at large wavevectors and should therefore be reliable 
for large widths where the inverse width provides a natural 
large wavevector cut-off. For widths w/£o > 2, the results are 
consistent with those of the MC simulations using compos- 
ite fermion trial wf's [p5[|. For 0.5 < w/£o < 2, the results 
are still comparable to those of the MC simulations, although 
they imply gaps which are rapidly increasing as w/^o ^ 0, in 
contrast to the results in Figure 1 1. 

Evidence, also based on CF trial wavefunctions, was pre- 
sented in which suggested a phase transition as a function 
of increasing width from incompressible states to compress- 
ible states at filling fractions v = /:'/(2/:i + 1). We have tested 
this theoretical prediction by computing the width dependent 
energy of the lowest energy excitation, which corresponds to 
the roton minimum. We have not analyzed the extrapolation 
to the infinite-system size limit for the roton minimum and 
present, instead, the variation with width of the roton min- 
imum energy for a system with fixed particle number. We 
show the results for v = 1/3, 2/5 and 3/7 in Figure |l^. We 
find that even for very large and unphysical widths up to hun- 
dreds of nanometers (corresponding to up to 50 £o), there is no 
evidence of the gap vanishing at any of these filling fractions. 
Instead we find that for such large width parameters the roton 
minimum scales simply as 1 /w. 



B. V = I state 



(0) , cos^A sin^* , 

where x = w/£o. The functional form ( ^^ is suggested by the 
following observations: The Haldane pseudopotentials V,„ for 

m > behave for w as ¥,„ w Vm^ + 0{w^) while for very 
large w they behave as Vm (log(w) /y/K + Um)/ w. Indeed, as 
will be seen below, the energy gaps decrease as 1 /w for very 
large values of the width w, implying that the logarithmic term 
cancels out in this limit. The values of the fitting parameters 
Eq \ (|), a and b are listed in Table pi for the filling fractions 
v = i 2 3 

^ 3' 5' T I 1 

The results presented in Figur e |1 1| are similar to those re- 
ported in p^]. The results of [[25[] were based on Monte 



Here, we present the results of calculations of the finite 
width effect on the energy gap of the mysterious V = 5/2 state. 
If the effects of Landau level mixing are neglected, it is suffi- 
cient to solve for the ground state of the electrons occupying 
the second - half-filled - Landau level and take the filled low- 
est Landau levels of spin-up and spin-down electrons as inert, 
i.e. unpolarizable. This problem is characterized by a filling 
factor v''' of the first excited Landau level of v''^ — 1 /2. It 
is customary to represent the system of electrons filling half 
the second Landau level by lowest Landau level wave func- 
tions but to take into account the interaction of electrons in the 
second Landau level by using the appropriate Haldane pseu- 
dopotentials of the second Landau level. Again, as for the 
computation of energy gaps at v = p/{2p + 1), there are es- 
sentially two ways to compute the energy. Either one may 
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Figure 12: The energy of the roton minimum as a function of width 
for systems with 11 particles at v = 1/3, 14 particles at 2/5 and 15 
particles at 3/7. For clarity, the energies of the roton are scaled by 
a factor 1/10 at v = 2/5 and 1/100 at v = 3/7. If there were to be 
phase transition to a compressible state the gap would have to vanish 
at some finite width. Instead we find clear evidence that, for large 
widths, the energy of the roton minimum (the lowest-lying excitation 
at fixed total flux) tends to zero as 1 /w with no suggestion of a phase 
transition. 



calculate neutral excitation (exciton) energies corresponding 
to a widely separated quasiparticle and quasihole pair, or one 
may calculate the energy of ground states containing a (frac- 
tionally) charged excitation. In the case of the v = 5/2 state, 
or equivalently at v'^^ = 1/2. there is the problem that elemen- 
tary charged excitations are predicted to occur only in pairs. 

The polarized ground state at v''^ — 1/2 occurs on the 
sphere when the number of flux units is 



2So = 2N-3, 



(27) 



and is thought to be described by a paired state, which may 
be of the Moore-Read pfaffian type [|l], H], However, 
great care is needed when analyzing excitation energies in 
these states on the sphere to avoid mistaking systems at con- 
ventional filling fractions Vp = p/{2p + I) or 1 —Vp for sys- 
tems at filling v''^ = 1/2. As we have discussed previously 
|[lC|], systems on the surface of a sphere exhibit degeneracies 
where, for a certain size A^, states with different filling factor 
coincide. This turns out to be a particularly severe problem 
in the sequence (p7[). Indeed, of the possible systems with 



N 
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Table III: Total flux in the ground (2So) and excited (25o ±1) states 
for systems at v'^' = 1/2 as a function of number of particles. A'. 
Where these states are aliased to conventional quantum Hall state 
ground states, we also show the corresponding filling fractions, Va- 
We note that the only sizes for which no aliases occur, are A' = 10,14 
and 18. Unaliased ground states occur in addition at A' = 8 and 16. 



up to 18 electrons, only five are not aliased with conventional 
fractional states, namely those with = 8, 10, 14, 16 and 18 
particles. Of these, the ones at = 8 and 16 have the problem 
that charged excitations of these states are aliased with ground 
states of conventional FQH states. Using these aliased states 
for a calculation of the energy gap at v^'^ = 1/2 would be mis- 
leading and would give rise to systematic errors. In Table [ll[ 
we list the relevant states and their aliases. We first show the 
energy of neutral excitations (exciton) with maximal angular 
momentum L^ax, corresponding to the largest possible sepa- 
ration of the quasiparticle and quasihole on the sphere. The 
angular momentum of this exciton is given by L,„ax — N/2 if 
N/2 is even, otherwise L,„ax — N/2 — I. In Figure the ex- 
citon energy for zero width, corrected for the Coulomb attrac- 
tion between quasiparticle and quasihole (A j /4 ( 1 /2 ) , equation 
(^), is plotted as a function of system size 1/A^ together with a 
linear fit in 1 /N to the data at A^ = 8, 10, 14, 16 and 18, cf. Ta- 
ble |l|. Like at V = 2/5, the exciton energy shows very large, 
and fluctuating finite size effects. Extrapolation to the bulk 
limit using a linear fit in 1 /N yields the result 



Av.-0.028-. 



(28) 



Alternatively, the energy gap can also be computed by cal- 
culating individually the energy of quasiparticle and quasihole 
excitations. The two quasiparticle state occurs at 250 — 1 while 
the two quasihole state occurs for 2Sq + 1 . Since in both cases 
the two excitations have the same charge {q — e/4 for the 
quasiparticle and q = —e/4 for the quasihole), one expects 
that the lowest energy state occurs when the two charges are 
maximally far apart, which demands maximum relative an- 
gular momentum, and consequently minimum total angular 
momentum on the sphere. Although one might have expected 
that this would imply L — Q for the ground state, as a result of 
symmetry, the angular momentum of the lowest energy states 
is L = N/2mod2, i.e. L = 1 for A^ = 10, 14, 18. The energy 
of these two-quasiparticle or two-quasihole states contains, in 
addition to the term A2^(v''') (equation^), the Coulomb in- 
teraction, AAg, of two quasiparticles separated by twice the 
radius R (the maximal separation on the sphere): 
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Figure 13: The exciton energy at v = 5/2 for zero width, cor- 
rected for the Coulomb attraction between quasiparticle and quasi- 
hole (Ai/4(l/2), equation is plotted as a function of system size 
1/N. 




(29) 



Combining the two contributions A2^(v''') and AA^ gives for 
the finite size coiTection term Q (A^) 
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(30) 



The gap calculation then proceeds by taking account explicitly 
of the finite size correction Cq{N) (equation as described 
for the cases at v = /?/ (2p + 1 ) in the previous section. 

In Figure |l^a), we show our results for the gap at v = 5 /2. 
In the top figure, half the sum of quasiparticle and quasihole 
excitation energies are plotted as a function of system size 
1 /N for different values of the width w. For zero width, re- 
sults for = 10, 14 and 18 are plotted, the sizes at which no 
aliasing effects occur. They lie almost exactly on a straight 
line in 1/A^. Extrapolation to the bulk limit yields 



A,/,«0.025-, 



(31) 



slightly lower, but consistent with the result (^8]) based on the 
exciton energies. Based on our previous experience with gap 
calculations at v = 1/3 and 2/5, we believe that also at v = 5/2 
the extrapolation based on individual quasiparticle and quasi- 
hole energies is more reliable. However, the exciton energy 
calculation certainly supports our conclusion that the quasi- 
particle and quasihole states at v = 5/2 contain two charged 
defects. Otherwise, there would be a factor of two difference 
between our extrapolated values A^^j (equation ^ and A5 /2 
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Figure 14: Energy gap at v = 5/2. The upper panel illustrates 
the size dependence of the energy gap for different values of the 
width parameter Q<wIIq< 1.17. The width parameters are vv/^q = 
0, 0.098, 0.195, 0.391, 0.586, 0.781, 1.17,with the topmost line re- 
ferring to the case w = 0, and the rest in sequence down to the low- 
ermost line with w/Iq = \.\1 . The extrapolations to the N limit 
assume that the slopes for the cases with w 7^ are the same as for 
the VI' = case. In the lower panel the gap values, extrapolated to the 
N ^ 0°, are plotted as a function of w/io. 



(equation^). Finally, in Figure [14|(b), the gap in the thermo- 
dynamic limit is plotted as a function of width w/Iq. These 
results indicate that the width effects reduce the gap at v = 5/2 
slightly. 

Very recently, Eisenstein et al [ p^ have investigated the 
V — 5/2 and 7/2 states in a sample of ultra-high mobility 
(/J 3.1 X lO^cm^/Vs). They determined an activation gap 
A™ 2 ~ 0.3 IK at V = 5/2 and A™ 2 ~ 0-07K at v = 7/2. Their 
sample had an electron density = 3 x 10^ Vcm^, which leads 
to a width w « 4nm. At v = 5/2, the field B = 4.96T cor- 
responds to a value w/£o « 0.35, while at v = 7/2 we get 
w/io « 0.30. According to the results shown in Figure |l^ 
The calculated gap values (see Figure for w/£o w 0.35 
and 0.30 are around 0.0220 and 0.0225?/e£o respectively. 
These lead to theoretical values for the gap of Ay 2 ~ 2.5K 
and Ay 2 « 2.1K. A disorder broadening of the order of 2K 
would explain the measured gaps of 0.3 IK and 0.07K. It is 
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Sample A (F = 


1.28±.13K) 


Sample B (F = 


2.1 ±.17 K) 


V 


B[T] 


A^' +r [K] 


A^ [K] 




B[T] 


A™ +F[K] 


At;[K] 


K 


1/3 


13.9 


9.03 


15.2 


0.59 


28.5 


13.2 


18.9 


0.70 


2/5 


11.6 


4.48 


6.9 


0.65 


23.8 


6.5 


8.6 


0.75 


3/7 


10.8 


3.23 


4.8 


0.67 


23.2 


4.5 


5.7 


0.79 


4/9 


10.4 


2.23 


3.6 


0.61 


21.4 


3.3 


4.9 


0.68 



V 


B[T] 


A?' [K] 


A^ [K] 


a;^' [K] 


1/3 


21.0 


10.5 


16.6 


13.5±.5 


2/3 


10.8 


6.5 


13.1 


10.7±.5 


5/3 


4.5 


1.0 


9.3 


7.9±.3 



Table IV: Comparison of the measured gaps, Ay , in samples A 
(nominal density 1.12 x lO" cm^^) and B (nominal density 2.3 x 
lO'^cm^^) reported in [ |l2| | with the gaps computed for a Coulomb 
interaction but taking account of the finite- width effects. Ay . We 
have added a constant field-independent shift, F, for each sample 
which we estimate by comparing the functional dependence of the 
gap energies as a function of filling fraction, v, predicted by CP the- 
ory with that found in experiment. The range quoted for F gives the 
maximum and minimum found when the constant C' in (MQ) varies 
between 4.11 (our estimate of C' for the pure Coulomb interaction) 
and 9. 



important to note that previous experimental values of the ex- 
citation gap at V = 5/2 have been much smaller ^8|]. For 
samples with density ris = 2.3 x lO'Vcm^ the gap at v = 5/2 
was 0.1 IK p^. In this case, the width is w w 4.5nm, and at 
the field B = 3.65T, we obtain w/Iq w 0.34, and a theoretical 
gap value of 2. IK. The factor of ~ 3 difference between the 
recently reported activation gap [^] and the earlier estimate 
| |4^ ] in samples with similar densities suggests that the activa- 
tion gap is affected strongly by sample quality, and is likely to 
be dominated by the effects of disorder. By comparing the gap 
at 5/2 to those at 7/3 and 8/3, and also at v = p/{2p + 1), Pan 
et al. also concluded that a disorder broadening of the order 
2K was to be expected. 



VII. EXPERIMENTAL GAPS 

Estimates of the gaps for fractional quantum Hall systems 
have been reported for GaAs heterojunctions [ pi] , [l2[ ] and 
more recently for metal-insulator-structures (MIS) using or- 
ganic (pentacene and tetracene) semiconductor layers [^]. 
The recent measurements on organic MIS structures are par- 
ticularly interesting given the different separation of energy 
scales to that found in GaAs. The dielectric constant in 
tetracene is in the range e «3 to 4 (compared to e wl2.7 in 
GaAs), the band mass is ^ 1.3mg (0.07 in GaAs) and the g- 
factor is close to 2 (0.44 in GaAs). The larger band mass and 
the smaller dielectric constant mean that, for samples with the 
same density, the ratio of interaction energies to the Landau- 
level splitting is much larger in the tetracene structures than in 
GaAs and hence that Landau level mixing effects are expected 
to be larger. The larger g-factor gives a larger Zeeman energy, 
and hence makes spin-reversed excitations less likely than in 
GaAs heterostructures. 

We have estimated the gaps atv = 1/3, 2/5 and 3/7 for the 
two samples A and B of JT^]. We take the quoted density of 
the samples and assume a depletion density — «s/5 (this 
is typical of these samples [|ri|, although the results are not 



Table V: The activation energies as deduced from the temperature 
dependence of the longitudinal resistivity at filling fractions v = 
1/3, 2/3, 5/3, A", reported in [ju]] are compared to our values for 
the gaps. Ay. For reference, we also show the calculated values of 
Willett et al. | [ll| ] in the last column. These authors fixed the deple- 
tion density n^^pi and hence the width parameter w by requiring that 
the solution of (llSl) correctly reproduced the experimentally mea- 
sured sub-band splitting. They estimated the finite width corrections 
on the basis of the model interaction (p2[) to give Ay . 



sensitive to the exact value of njepi)- /,From the results in Fig 
we estimate the standard deviation of the density distribu- 
tion and this leads directly to an estimate of the gaps (see Fig 
[TTI ). We compare our results with those of the two samples 
A and B of |^4| in Table [V The effects of impurity scat- 



tering have been taken into account empirically by assuming 
that the levels are broadened by a field-independent broaden- 
ing F. This assumption has not been theoretically justified. 
However, for the purpose of comparison we have reanalyzed 
the results of [ p^ under this assumption by fitting them to the 
functional form predicted by CF theory, i.e. including the log- 
arithmic corrections (see |l^ to extract the broadening F. We 
find that the gaps measured are consistently between 60 and 
70% of what we predict after taking account of finite thick- 
ness effects. This is consistent with the results of |M, Ga], 



coiTecting the error of Reference 1 23 1. 

The results reported in [ pl| relate to filling fractions p/3, 
where p = 1, 2, 4 and 5 and were interpreted on the assump- 
tion that the ground states and gaps were all maximally spin- 
polarized states within the lowest Landau level for a sam- 
ple with density — 1.65 x lO" cm^^ and mobility 5 x 10^ 
cm^/V sec (to be compared with 6.8 and 12 x 10*' cm^/V sec 
in samples A and B in [p^). The authors of [|ll]] solved ( p3| ) 
numerically for the sub-band wf, ^(z), choosing the depletion 
density njgpi to reproduce the experimentally observed sub- 
band splitting. As as a result we have a more precise estimate 
of width of the wf in the lowest sub-band than we have been 

We have converted their 



able to make for the samples of 1 12|. 
estimate of the width to a standard deviation w and estimated 
the gaps at the relevant filling fractions. The results are pre- 
sented in Table M. We note that the measured values A'" of 
the gap at V = 1 /3 and 2/3 are both larger than our theoretical 
values by about the same amount F w 6K. This might serve as 
an estimate of the broadening. The authors of [O] also esti- 
mated the gap reduction on account of finite thickness effects 
based on the exact diagonalizations of six particle systems re- 
ported in (32 ) and we include these estimates Ay in the Table. 
The reduction of the gaps found in the earlier finite-size stud- 
ies was significantly larger than what we obtain (Section 3). 
It may have resulted from estimating the gap reduction using 
systems which were too small, or inaccurate extrapolation to 
the thermodynamic limit. 
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It is clear from both Tables IV and ^ that the discrepancy 
between measured gaps and computed gaps is significant. 
This discrepancy may be due to Landau level mixing, spin- 
reversed excitations and to impurity effects not accounted for 
by the use of the field-independent broadening F used in Table 



IV) Estimates in [gjj based on diagonalizations of up to only 
5 particles in a torus geometry suggested reductions of the gap 
(identified with the zone boundary exciton) as a result of Lan- 
dau level mixing of between 12% and 17% were possible at 
V = 1/3 in a magnetic field at lOT for a pure Coulomb inter- 



action. These should scale as {e /tto) /H&c ^ 1/vB- On this 
basis the reduction at a field of 28. 5T for sample B at v = 1/3 
would be at most 10%. However, as the matrix elements be- 
tween Landau levels of the effective interaction V{r\ — ^2!), 
which diverges only logarithmically as r ^ 0, will be signifi- 
cantly smaller than those of the bare Coulomb interaction, the 
reduction of the gap due to Landau level mixing in these sam- 
ples should, in fact, be significantly smaller than this figure of 
10% and is probably negligible. Clearly, a new study along 
the lines of [pQ] (which actually looked at the harder problem 
of Landau-level mixing at v = 5/2 for systems with a partially 
filled second Landau level), taking account of the finite-width 
of the subband wf, would make for significantly more accurate 
estimates of Landau level mixing effects. 

We should also consider the role of excitations involving 
spin reversals. The gap at v = 1/3, corresponding to the cre- 
ation of a quasihole with no spin-reversal and a quasiparti- 
cle with one spin-reversal, was estimated in [|[ |T]| using ex- 
tensive Monte Carlo simulations of trial wf 's. Gap estimates 
for various combinations of quasihole and quasiparticles com- 
bined with spin-reversals based on exact diagonalizations of 
small systems (up to six particles) were reported in [p9||. For 
the case of the pure Coulomb interaction and ignoring the Zee- 
man energy the gap to create spin-reversed excitations was 
around 60% smaller than the spin-polarized gaps for systems 
at V — 1/3. When the Zeeman energy (in GaAs) is taken into 
account this suggests that, for a pure Coulomb interaction, the 
spin-reversed excitation would have a lower energy for sys- 
tems at V = 1 /3 if the magnetic field were smaller than ^ 7T. 
This is above the fields at which the 4/3 and 5/3 states were 
observed in JTl] ] (see Table and may account for the larger 
discrepancy seen at these filling fractions than at v = 1/3 or 
2/3. 

The estimate of 7T, as the field below which the spin- 
reversed excitation drops below the spin polarized excitation, 
is well below the fields in Table |^ making it unlikely that 
spin-reversed excitations are involved at these filling frac- 
tions. Although the explicit estimate of the spin-reversed ex- 
citation was made for a system at v = 1/3 it is unlikely that 
the discrepancy at other filling fractions will be larger. This 
is because the difference between the spin-polarized and spin- 
reversed quasiparticle energies should be largest at v = 1/3, 
where it is possible to construct a spin-reversed quasiparticle 
state which is a zero-energy eigenstate of the hard-core po- 
tential (with only the pseudopotentials Vo and Vi non-zero). 
For the case of the Coulomb interaction, its energy is con- 
trolled by the size of the pseudopotential V2, while the energy 
of the spin-polarized quasiparticle is determined by the larger 



Vi . However, as V2 is reduced less by finite width effects than 
Vi (see Figure the spin-polarized quasiparticle will be sta- 
bilized with respect to the spin-reversed excitation by finite- 
width effects 1^]. 

The results for the activation gaps at v = 1 /3 and 2/5 in lay- 
ers of tetracene reported in | p9| ] can also be compared with our 
numerical results. By simultaneously varying the gate voltage 
and magnetic field the gaps could be tracked as a function of 
the ratio w/Iq for a range of fields 3 < < 9. One in- 
triguing feature of these organic layers is that the ratio of the 
Coulomb interaction to Landau level spacing, 
is approximately 30-40 times larger than in the GaAs samples 
for systems at the same magnetic fields. 

We have computed the width, w, of the subband wavefunc- 
tions in the tetracene samples of | p9| ] excluding the effects of 
image charges using the zG trial wavefunction and found that 
w varies between 11 A at a density iis = 0.1 x lO'' cm^^ and 
lA for iis — Sy. 10"cm^^. The effects of the image potential 
will be to reduce the width still further At all the densities, at 
which the gaps were measured, w/^o < 0.1 and so the effects 
of the finite width of the wavefunction on the gaps in these 
samples are small (see Fig pTj) a nd significantly smaller than 
were reported in Reference However, our calcula- 

tions summarized in Fig |ri| are a more accurate reflection of 
finite width effects than the old formula of [ pO| ] used in [^] . 
We also note that the widths we obtain are about half the order 
of magnitude {w ^ 35 A) quoted in [p^. 

In order to compare the results of our calculations with the 
results of the experiments on the tetracene MIS structures, we 
need to take account of the large difference between the di- 
electric constant of the alumina insulating layer (e ~ 9.8) and 
the value for tetracene (e ^ 3.5). Nearly all the charge density 
is within 3w of the interface. This is significantly less than the 
average interparticle spacing or ion disk radius (a — yjlfwlo), 
which for these samples varies from \5QA to 360A depending 
on density. We have used the simplest approximation which 
treats the 2D electron gas as if it were trapped at the inter- 
face of the alumina and the semiconductor The effective di- 
electric constant is then just the average for the two materials. 
(Corrections to this, taking account of the actual displacement 
of the charge away from the interface, would involve image 
charge effects and give rise to a change in the functional form 
of the effective interaction between particles [^.) In Figure 
[T^ we compare the results we obtain with the measured val- 
ues reported in [ p^ ] . We show calculated gaps as a function of 
magnetic field ignoring the finite width of the charge distribu- 
tion. 

The difference between the computed and measured values 
of the gaps in Figure |l5| at v = 1/3 and 2/5 are remarkably 
small. Although there is some uncertainty associated with 
the computed gaps arising from the simple treatment of the 
large difference in dielectric constant of tetracene and alu- 
mina, there is surprisingly little evidence of large Landau level 
mixing (LLM) or disorder-related effects at these filling frac- 
tions. The ratio of the Coulomb energy scale e^/e^o to the 
cyclotron energy in tetracene is ^ 93 /\/B or ~ 42 at B = 5T 
and LLM should be significant and might even be expected to 
be dramatic. When the ratio of these two energy scales is this 
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Figure 15: Gap energies at v = 1/3 (upper panel, unfilled sym- 
bols) and 2/3 (upper panel, filled symbols) and 2/5 (lower panel) 
in tetracene samples from |^] compared with numerical estimates. 
Triangles and squares denote experimental results for different sam- 
ples. The solid line in each panel shows the gap given by the quoted 
formula which is valid in the zero thickness limit (as discussed in the 
text the effects of the non-zero width are negligible in these samples). 
The effective dielectric constant is taken as the average e = 6.65 of 
reported values for tetracene £ = 3.5 and the value e = 9.8 for the 
alumina insulating layer. The dashed line in the upper panel is an 
estimate of the energy gap for spin reversed excitations (see text). 



large a perturbative treatment of LLM effects may not even be 
possible. (We have assumed the same effective dielectric con- 
stant as used in Fig 15 ). Even though the mobilities are not as 



high in the tetracene MIS structures as in the GaAs-GaAlAs 
heterostructures |l^, the agreement between calculated 
and measured gaps suggest that there are not any strong ef- 
fects of disorder scattering either. 

For the systems at v = 2/3 (filled symbols in upper panel 
of Fig |5]), it was suggested in [ p^ that the change in slope at 
around 6.5T was related to a transition from a polarized state 
at high fields to a state which was not fully polarized at lower 
fields. This seems unlikely. At a transition with a disconti- 
nuity in polarization (first order transition), there would nor- 
mally be a discontinuity in the gap rather than a discontinuity 
in its gradient with magnetic field, see for example | ]54j ] . There 
has been one report of a transition from a polarized to partially 
polarized state at v = 2/3 in GaAs heterostructures without 
any discontinuity of the gap [|5^]. However the corresponding 
transition would be expected to occur in the tetracene samples 
at around 1 .7T well below the range of fields of Fig |l5[ Even 
if there were no (or only a small) discontinuity in the gap. 



the change in slope would normally be in the opposite sense 
to the one reported (see Figure |l5|). The lowest-lying excita- 
tions from a partially polarized (or unpolarized) state would 
be expected to involve spin reversals which increased the to- 
tal aligned spin (rather than reduced it) and thereby gained a 
reduction in Zeeman energy. On the other hand, excitations 
from the fully polarized state, either decrease the total spin or 
leave it constant. As a result there would either be a contribu- 
tion to the energy of the excitation from the Zeeman energy, 
which was positive and increasing as a function of field, or 
no contribution. In either case, the gap would be expected to 
grow faster with field in the high field (fully polarized) state 
than in the low field state but not not more slowly as reported 
in [p^]. This is what was observed for the transition seen at 
v-8/5 @. 

An alternative explanation of the results at v = 2/3, as- 
sumes a fully polarized ground state and identifies the change 
in slope at B « 6.5T with a change in the nature of the low- 
est lying excitations. For B > 6.5T, the lowest energy exci- 
tations would be within the fully spin-polarized sector, while 
for B < 6.5T they would involve a spin-reversal. We can make 
a rough estimate of the energy of the spin-reversed excitation 
gap ignoring Landau level mixing (LLM) as follows. Previous 
estimates of the Coulomb energy of a spin-reversed quasipar- 
ticle energy put it at around 55% of the energy of the spin- 
polarized quasiparticle ||^]. The Coulomb energy of a spin- 
reversed quasihole is unlikely to be much lower than that of 
the spin-polarized hole, while the additional Zeeman energy 
will make this excitation unfavorable. We therefore take for 
the value of the quasihole energy that of the spin-polarized 
hole. The results for this 'spin-reversed' energy gap are shown 
as a dashed line in Figure ( [l5| ) and are seen to be quite close 
to the observed data points, although the difference between 
our results for the 'polarized' gap and the spin reversed' gap 
is small. Our estimate of the spin-reversed excitation applies 
both at V = 1 /3 and v = 2/3 as we have neglected LLM, which 
allows electrons to make virtual transitions to other Landau 
levels and thereby screen the interaction in the lowest Landau 
level. These effects would be larger for systems at v = 2/3 
than at V = 1 /3 in the same magnetic field and could ex- 
plain why the the spin-reversed excitation lies below the fully 
polarized excitation up to higher magnetic fields at v = 2/3 
than at V = 1/3. A tilted field experiment [|4| would be one 
method to determine whether our identification of the change 
of slope in the gap with field with a change in polarization of 
the lowest-lying excitation is correct. 

The apparent absence of a significant reduction of the gap 
in the tetracene MIS structures (// < 2.5 x 10^ cm^A'^s) on ac- 
count of disorder, given its importance in the ultra-high mobil- 
ity (12.8 X 10^ cm^A's) GaAs heterostructures, is puzzling. It 
suggests that the activated gap measured in transport measure- 
ments is affected by disorder in different ways in heterostruc- 
tures and MIS structures. In the heterostructures, the disorder 
scattering is that of the ionized silicon donors which were in 
a layer about 8004 from the electi-ons |Q. In the MIS struc- 
tures, on the other hand, the doping is controlled by a capaci- 
tance (~ 1 30 nFcm^^ p^]) with the backgate of order microns 
from the carriers. Here the disorder scattering is likely to be 
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that of neutral defects. It is possible that, in the heterostruc- 
tures, the activation studies do not measure directly the energy 
to create a quasiparticle quasihole pair from the ground state, 
but rather the energy to excite quasiparticle (or quasiholes) 
out of bound states in the potential of the (charged) impurity 
distribution. 



Vm. CONCLUSIONS 

We have used diagonalizations of the Hamiltonians for fi- 
nite size systems on a sphere to obtain estimates of the gaps 
at filling fractions in the Jain sequence v = 1/3,2/5,3/7 and 
4/9 and at v = 5/2. We have emphasized how taking account 
properly of the systematic contributions to the excitation en- 
ergy from the charge redistribution on the sphere in excited 
states is essential if one is to obtain accurate estimates of the 
gaps in the thermodynamic limit. Our results for the gaps are 
smaller than earlier estimates based on finite-size studies (for 
V = 2/5 and 3/7 0]) and those based on the study of trial 
wavefunctions (for v = 2/5, 3/7 and 4/9 [|l9ll). This differ- 
ence is important, as estimates of the gap as a function of v 
provide the most direct numerical estimates the effective mass 
of CP's 1 15, O]. Our new results are consistent with the CF 



picture provided the logarithmic corrections to the effective 
mass are taken into account and are not well described by as- 
suming a filling factor independent effective mass (see Fig^). 

We have shown that Gaussian and the zx Gaussian (zG) 
variational functions accurately describe subband wavefunc- 
tions and yield subband energies and lowest Landau level 
pseudopotentials, which are essentially indistinguishable from 
those obtained by solving for the subband wavefunctions ex- 
actly by direct numerical integration. These trial wavefunc- 
tions offer a significant improvement over the standard Fang- 
Howard (FH) form, which overestimates the standard devi- 
ation of the charge distribution, w, by as much as 50% de- 
pending on electron density, n^, (see Fig 0). The lowest Lan- 
dau level pseudopotentials, which are the starting point for the 
study of the fractional quantum Hall gaps, turn out to be accu- 
rately determined using any of the three trial forms (Gaussian, 
zG or FH) once w is known (see Fig This offers a huge 
simplification over the previous ab initio approaches which 
used numerical integration to find the subband wavefunction 
and tables of pseudopotential parameters p^. We have also 
computed the variation of the gaps at fractionally quantized 
Hall states as a function of width of the subband charge distri- 
bution. The results are parametrized in Eq |2^ and Table ^ 

We have compared our computed gaps with measured acti- 
vated gaps. We have found that, even after taking account of 
disorder broadening of states, the measured activation gaps in 
GaAs heterostructures are only around 60% of the computed 
gaps for the filling fractions v = 1/3,2/5,3/7 and 4/9 (see 
Table |^). This is to be contrasted with the activated gaps at 
V = 1/3 and 2/5 reported in tetracene MIS structures, which 
turn out to be reasonable agreement with computed gaps (see 
Fig |T^). We have suggested that the relationship between the 
computed gap and the measured activated gap may be dif- 
ferent depending on the type of disorder in the samples. In 



the GaAs heterostructures the charged donor ions, which are 
the main scattering centers, are only around 800A from the 
quantum Hall layer, and this could lead to local variations 
in the energy required to excite quasiparticle-quasihole pairs 
with the lowest excitation energies controlling the activated 
transport. On the other hand, the backgate in the MIS struc- 
tures is of order microns from the quantum Hall layer, and the 
main scattering centers are likely to be neutral. These are less 
likely to affect the energy to excite quasiparticle-quasihole 
pairs and the gap controlling activated transport should then 
be close to the true thermodynamic gap as we have found. 
For filling fractions v = 2/3 and 5/2, the reported activated 
transport gaps in GaAs heterostructures are only around 10% 
and 5% respectively of the values we compute, although we 
have not attempted to account for disorder broadening in these 
cases. However, the experimental evidence suggests a disor- 
der broadening which is comparable at v = 5/2 with the com- 
puted gap [Est] so the large discrepancy is to be expected. 
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Appendix: ZdS Interaction 

One previous attempt to model finite width effects used the 
'model interaction' [ |20| | 



Vzdsir) 



72' 



(32) 



which introduces a width parameter t. We have found that this 
model interaction cannot reproduce accurately the variation 
with L of the Haldane pseudopotential parameters for a sample 
with finite width with the errors significantly increasing as the 
width increases. Figure |l^. 

The reason for this is probably the unphysical nature of this 
interaction as a model for electrons in a heterostructure or 
quantum well interacting via the Coulomb interaction. Taking 
the Fourier transform of (|l4|) and using the convolution theo- 
rem, one can show that it is not possible to construct a densit 



nsity 

4i 



distribution, n for which the effective interaction (sei 
is Vzdsi'')- This is essentially because Vzds{>') is the Coulomb 
interaction of two (distinguishable) particles confined to sepa- 
rate planes a distance t apart and, as such, misses the In r found 
for small r and large widths for all realistic density distribu- 
tions. However, many of the results obtained on the 
basis of the effective interaction are still valid if interpreted 
carefully. 

For this purpose, we incorporate an overall scaUng factor of 
the interaction A^, 



vzdsir) ^N- 



1 



(33) 
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Figure 16: Errors in the Haldane pseudopotentials Vl computed us- 
ing the Zhang-DasSarma model interaction (^2|). The comparison is 
with the values reported in [p^]. The variational parameters are de- 
termined such that V] is correctly reproduced. Note in particular the 
large errors for Vq and the slow decay of the error with increasing L 
for the low-density sample. 



with N = I giving the original interaction (|32|). The gap ener- 
gies and relative stability of fractional quantum Hall states in 
the principal Jain sequence are determined principally by the 
first two Haldane pseudopotentials for odd angular momen- 
tum Vi and V3 . In figure O we show the values of t required in 
^ to match the values of Vi and of V3 to those obtained using 
the variational Gaussian wf as a function of the width param- 
eter assuming = 1 . It is clear that it is not possible to find 
a value of t which gives both V\ and V3 correctly. If we allow 
and t to vary then both Vi and V3 can be correctly repro- 
duced by the effective interaction in (^3|). The results are also 
shown in Fig |l7[ Changing means that the asymptotic be- 
havior of the pseudopotentials at large angular momentum is 
not reproduced correctly. However, as the gaps and stability of 
the incompressible states in the Jain sequence are determined 
principally by the pseudopotentials Vi and V3 this should not 
be a problem. If phase transitions between spin-singlet and 
polarized states (e.g. at v = 2/5) are of interest, it is obviously 
possible to correctly represent the in this case most important 
pseudopotentials Vi and V2 for angular momenta L = 1 and 2, 
by appropriate choice of and f . 
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Figure 17: The width parameter t (lower three curves) in the ZDS 
model interaction in(p3|) and the normalization parameter A' (top 
curve), as a function of the Gaussian width w of the density distri- 
bution, ^(z)^. For the two curves labelled A' = 1, the normalization 
is held equal to one and the parameter t is chosen so that the Hal- 
dane pseudopotential Vi (lower curve) and V3 (upper curve) for the 
interaction in ^ are equal to the values obtained from ( p^ using the 
Gaussian variational wf 's. The dashed-dotted curves show the values 
of t and A' required to reproduce both Vi and V3 /V\ correctly. 
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